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ABSTRACT 



Aims. Hyper-Luminous Infrared Galaxies (HLIRGs) are the most luminous persistent objects in the Universe. They exhibit extremely 
high star formation rates, and most of them seem to harbour an Active Galactic Nucleus (AGN). They are unique laboratories to 
investigate utmost star formation, and its connection to super-massive black hole growth. X-ray studies of HLIRGs have the potential 
to unravel the AGN contribution to the bolometric output from these bright objects. 

Methods. We have selected a sample of 14 HLIRGs observed by XMM-Newton (type 1, type 2 AGN and starburst), 5 of which 
are candidates to be Compton-thick objects. This is the first time that a systematic study of this type of objects is carried out in the 
X-ray spectral band. Their X-ray spectral properties have been correlated with their infrared luminosities, estimated by IRAS, ISO 
and sub-millimeter observations. 

Results. The X-ray spectra of HLIRGs present heterogeneous properties. All our X-ray detected HLIRGs (10) have AGN-dominated 
X-ray spectra. The hard X-ray luminosity of 8 of them is consistent with a pure AGN contribution, while in the remaining 2 sources 
both an AGN and a starburst seem to contribute to the overall emission. We found soft excess emission in 5 sources. In one of them it 
is consistent with a pure starburst origin, while in the other 4 sources it is consistent with an AGN origin. The observed X-ray emission 
is systematically below the one expected for a standard local QSO of the same IR luminosity, suggesting the possible presence of 
absorption in type 2 objects and/or a departure from a standard spectral energy distribution of QSO. The X-ray-to-IR-luminosity ratio 
is constant with redshift, indicating similar evolutions for the AGN and starburst component, and that their respective power sources 
could be physically related. 
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1. Introduction 

Ultra-luminous Infrared Galaxies (ULIRGs) are a class of galax- 
ies with bolometric luminosity L/s > 10 '^L p, dominated by the 
emis sion in the infrared (IR) waveband (see lSanders & Mirabell 
11996 for a complete review). They are, together with opti- 
cal quasars, the most luminous objects in the Local Univ erse. 
ULIRGs are rare in the Local Universe dSoifer et al.lll987h . but 
large numbers are detected instead in deep-IR surveys, and are 
a fundamental constituent of the high redshift galaxy population 
dSmail et al.lll997l: iGenzel & Ces arskv 2000; Franc eschini et al.l 
1200 lb . They are powered by Active Galactic Nucleus (AGN) 
and/or st arburst (SB) trigger ed by mergers of gas-rich spiral 
galaxies dVeilleux et al.ll2002h . Optical spectroscopic studies re- 
veal that the fraction of ULIRG s hosting an AGN i n crease s 
with increasing IR luminosity (IVeilleux et al.l 1 19951 Il999l) . 
Furthermore, the fraction of Seyfert 1 to Seyfert 2 ULIRGs in- 
creases with IR luminosity. 

It has been proposed that ULIRGs at high redshift 
could be the origin of massive elliptic a l and SO galaxies 
dFraiiceschini et al.lll994l; iLillv et al.lll999l: iGenzel & Cesarskv! 
I200d l). An important fraction of stars in present day galaxies 
would have been formed during these evolutionary phases. 

Observations in the X-ray band are a powerful tool to 
disentangle the AGN contribution to the bolometric luminos- 
ity from the ULIRGs. X-rays studies of ULIRGs have con- 
firmed their composite nature (powered by AGN and/or star- 
burst), with indications for a predominance of the SB over 
the AGN phenomenon, even when observed in hard X-rays 
dPranceschini et al.ll2003"l) . 



The brightest end of the ULIRG distribution is occupied 
by the Hyper-Luminous Infrared Galaxies (HLIRGs, L/r > 
IO'^Lq). They are among the most luminous objects in the 
Universe, although the origin of this luminosity is still uncertain. 
This population exhibits extremely h igh star formation rates, and 
most seem to also harbour an AGN dRowan-Robinsoij2000l) . 

The source and trigger of the emission from HLIRGs have 
been discussed since their discovery. Initially the extreme lumi- 
nosity of HLIRGs was attributed to gravitational magnification, 
but Hubble Space Telescope (HST) observations discovered that 
only a minority of these galaxi es (~15-20 per cent ) have been 
misclassified owing to lensing ( Farrah et al.ll2002"bl) . CuiTently, 
there are three main hypotheses on the nature of these objects: 

a) They could be simply the high luminosity tail of the ULIRG 
population, and therefore their power sources are proba- 
bly triggered by galaxy mergers. Though HST has revealed 
some merging systems among HLIR Gs, there is a signific ant 
fraction of them in isolated systems dFarrah et al.ll2002bl) . 

b) Assuming that the majority of the rest-frame far infrared 
(FIR) and sub -milUmeter (sub-mm) emission comes from 
star formation dVerma et al.ll2002l; iFarrah et al.ll2002a[) . their 
estimated star formation rates (SFR) are > lOOOMoyr \ the 
highest for any object in the Universe. HLIRGs could be very 
young g alaxies going through th eir major episode of star for- 
mation (lRowan-Robinsonl2000l) . 

c) They may be a completely new class of objects, where 
the IR emission is originated via some other mechanism: 
e.g., a transient IR -Iuminous phase in quasar evolution 
dFarrah et al.ll2002bl) . 
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Table 1. Hyper-luminous galaxies observed by XMM-Newton with redshift less than ~ 2. 



Source 


Type" 


IR-Model" 


RA 


DEC 


z 








AGN / SB 










[erg s"' cm"^] 


IRAS 00182-7112 


QSO 2 (CT) 


0.35 / 0.65 


RROO 


00 20 34.7 


-70 55 27 


0.327 


46.49 


IRAS F00235+1024 


NL SB (CT) 


0.5/0.5 


F02a 


00 26 06.5 


+ 1041 32 


0.575 


<46.76 


IRAS 07380-2342 


NL 


0.6 / 0.4 


F02a 


07 40 09.8 


-23 49 58 


0.292 


46.56 


IRAS 09104+4109^ 


QSO 2 (CT), cD 


1/0 


RROO 


09 13 45.4 


+40 56 28 


0.442 


<46.42 


PG 1206+459 


QSO 


1/0 


RROO 


12 08 58.0 


+45 40 36 


1.158 


47.20 


PG 1247+267 


QSO 


1/0 


RROO 


12 50 05.7 


+26 31 08 


2.038 


47.70 


IRAS F12509+3 122 


QSO 


0.6 / 0.4 


F02a 


12 53 17.6 


+31 05 50 


0.780 


<46.86 


IRAS 12514+1027 


Sy2 (CT) 


0.4 / 0.6 


RROO 


12 54 00.8 


+ 10 11 12 


0.3 


46.18 


IRAS 13279+3401 


QSO"* 


0.7/0.3 


F02a 


13 30 15.3 


+33 46 29 


0.36 


46.58 


IRAS 14026+4341 


QSO 1.5'' 


0.6 / 0.4 


F02a 


14 04 38.8 


+43 27 07 


0.323 


46.26 


IRAS F14218+3845 


QSO 


0.2 / 0.8 


F02a 


14 23 55.0 


+38 32 14 


1.21 


47.80 


IRAS F15307+3252 


QSO 2 (CT) 


0.7/0.3 


V02 


15 32 44.0 


+32 42 47 


0.926 


<47.07 


IRAS 16347+7037 


QSO 


0.8/0.2 


F02b 


16 34 28.9 


+70 31 33 


1.334 


47.42 


IRAS 18216+6418" 


QSO 1.2, cD 


0.6 / 0.4 


F02a 


18 21 57.3 


+64 20 36 


0.297 


46.49 



" NL: narrow-line objects; Sy2: S eyfert 2. C ompton-thick candidates are labeled as CT. Spectral classification from Rowan-Robinson' ('2OO0j), 
except IRAS F00235+1024 ( Ver ma etal][2002 ), IRAS 14026+4321 ( Wan g et al. 20063 and IRAS 18216+6 418 ( Ver on-Cetty & Ve ron 200d). 
IRAS 00182-7112 has been classify as type 2 source using the optical emission lines from lArmus et al. and the diagnostic diagrams from 
lOsterbrockl {19891 chap. 12). 

* Fraction of the IR e mission originated in AG N and/or SB. See Sect. |2]for details. Data from: (RROO) iRowan-Robinsoi] |2000|. (F02a) 
iFarrah et aU2002al . ( V02) I Verma etaUl2002l . (F02b) lFm^irn]|2002bl . 
Source in cluster. 

Not detected in X-rays. We use the optical data to classify this source as QSO. See Sect.|2]for details. 



The state-of-the-art X-ray, IR and sub-mm observations sug- 
gest that HLIRGs are powered by dus t-enshrouded black hole 
accretion and viol e nt star form ation ('Rowan- Robinson 2000; 
Verma et alJ 120021: iFarrah et ah 2002a; Wilman et al, 2003; 



Iwasawa et alj 2005h . pointing to a mixture of the (a) and (b) 
hypotheses. Recent observations of individual HLIRGs with 
Chandra and XMM-Newton show that the IR emission of these 
objects could be powered by buried quasars through dust re- 
radiation. The nuclear sour ce is heavily obscured, reaching the 
Compton-thick (CT) limit (llwasawa et al.ll200U IWilman et aP 
I2OO3I: I Iwasawa et a n I2005h . It has also been suggested that 
galaxy mergers in an over-density region may be a nec- 
essary condition for the formation of this class of sources 
(llwasawa et al. I I2005I) . However this suggestion is in contradic- 
tion with the signifi cant fraction of isolated HLIRGs found by 
IFarrah et alJ(l2()02b ^. 

AGN research has taken a central role in the study of 
galaxy formation and evolution since the discovery that most lo- 
cal spheroidal galaxy components (elliptically galaxies and the 
bulge of spiral galaxies) contain a super-massive black hole with 
the mass being proportional to the stellar mass of the galaxy 
spheroid (Magorrian et al. 1998; McLure & Dunlop 2002), 

These facts are most easily explained if the formation of the 
super-massive black hole and the spheroid of galaxies were co- 
eval, i.e., the black hole was built up by accreti on of the same gas 
that rapidly formed th e stars of the spheroid dPage et alj|2004l; 
iDi Matteo et"aLll2005h . 

Additional evidence for the coeval hypothesis is the simi- 
larity between the evolution of cosmic star formation and lu- 
minous AGN activity, since both were much higher in the past 
up to redshift ~ 2, than in t he present day dFranceschini et alJ 
1 19991: ISilvermanetani2005l) . As HLIRGs could represent the 
most vigorous stage of galaxy formation, they are unique labo- 
ratories to investigate extremely high stellar formation, and its 
connection to super-massive black hole growth. 

The main objective of this paper is to determine the relative 
contribution of SB and AGN emission to the bolometric lumi- 



nosity in HLIRGs, and the interplay between those two phe- 
nomena, as well as how their relative contribution varies with 
cosmic time. In our study we have also compared the properties 
of HLIRGs to those of ULIRGs. 

We have used X-ray observations of HLIRGs from XMM- 
Newton and also public data from the Chandra archive. In 
Sect. |2| the sample selection is explained. Data reduction and 
spectral analysis of each source are described in Sect. [3] Results 
are discussed in Sect.|4| Section |5] summarizes our conclusions. 

The Wilkinson Microwave Anisotropy Probe {WMAP) con- 
cordance cosmology has been adopted in the paper: Hq - 70 km 
s ' Mpc ', 0,„ = 0.27, Qa = 0.73 (.Spergel et al...2003.) . 

2. The HLIRG sample 

Our sample of HLI RGs has been selected from the 
Ro wan-RobinsonI (l2000i) sample of 45 HLIRGs. The sources in 
this mother sample can be classified in four sub-samples: (1) 
objects found from direct optical follow-up of 60 /im or 850 
surveys; (2) sources found from comparison of known quasar 
and radio galaxy lists with 60 catalogues, or using IR color 
selection (biased to AGN); (3) sources selected from sub-mm 
observations of known high-redshift AGN; and (4) known 
luminous IR galaxies not included in the previous subsamples, 
but satisfying L/s > IO'-'Lq. The first one is a flux limited 
sample, unbiased towards AGN; the sources in the second and 
third sub-samples have been selected in order to host an AGN, 
and therefore suffer from selection effects. 

We have chosen all the HLIRGs with public data in the 
XMM-Newton Science Archive (XSA), as of December 2004. 
We included also observations of five sources from OBS-ID 
030536 by our group (see Table O. Then we constrained the 
resulting sample to those sources with redshift less than ~2, to 
prevent a strong bias due to the presence of high-z QSOs. In our 
final sample there are 8 sources which are included in the first 
Ro wan-RobinsonI (12000) subsample, 4 sources are in the second 
subsample, 1 source is in the third and 1 is in the fourth one. 
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Source Net. exp. time [ks] Ext.radius [arcsec] Source counts" Obs. date Filter 





PN 


MOSl 


MOS2 


PN 


MOSl 


M0S2 


PN 


MOSl 


MOS2 






IRAS 00182-7112 


9.1 


- 


- 


19 


- 


- 


134 ± 15 


- 






2003-04-17 


Thin 


IRAS F00235+1024 


14.4 


_ 


_ 


35 


_ 


_ 


< 30 


_ 






2001-01-10 


Thin 


IRAS 07380-2342 


4.7 


11.4 


11.6 


30 


40 


40 


<45 


- 






2005-10-13'- 


Medium 


IRAS 09104+4109* 


9.2 


12.6 


12.5 


20 


20 


20 


5544 + 75 


2245 ± 47 


2332 


±48 


2003-04-27 


Medium 


PG 1206+459 


5.9 


6.9 


6.9 


26 


32 


29 


731 ±32 


188 ± 16 


193 : 


fc 16 


2002-05-11 


Thin 


PG 1247+267 


19.4 


25.5 


26.6 


35 


35 


35 


5132 ±74 


1646 ± 42 


1694 


±43 


2003-06-18 


Medium 


IRAS F 12509+3 122 


11.9 


14.6 


14.1 


25 


25 


25 


508 ± 24 


156 ± 13 


139: 


fc 13 


2005-12-1 r 


Thin 


IRAS 12514+1027 


16.7 


20.0 


18.2 


22 


40 


40 


105 ± 22 


28 ± 13 


7± 


12 


2001-12-28 


Thin 


IRAS 13279+3401 


24.6 






25 






< 36 








2006-01-18'' 


Medium 


IRAS 14026+4341 




6.6 


5.6 




35 


35 




< 30 


< 27 


2005-11-26'- 


Thin 


IRAS F14218+3845 


11.5 


15.2 


15.1 


30 


34 


31 


550 ± 29 


185 ± 17 


176: 


fc 16 


2003-08-01 


Medium 




2.3 


7.3 


7.0 


30 


30 


30 


91 ± 14 


109 ± 20 


79 ± 


18 


2005-06-07'- 


Medium 


IRAS F15307+3252 


9.4 


11.5 


12.0 


22 


40 


40 


97 ±21 


30± 12 


42 ± 


13 


2002-07-30 


Medium 


IRAS 16347+7037 


12.9 


15.6 


15.8 


30 


56 


48 


11172 ± 106 


3546 ± 67 


3785 


±66 


2002-11-23 


Medium 


IRAS 18216+6418* 


0.5 


2.9 


3.3 


20 


20 


20 




7402 ± 86 


8481 


±92 


2002-10-16 


Thin 



" Total source counts in the 0.2-10 keV band. 

* Chandra data are also used in the study of this source. See Section [331 for details. 
' Data from OBS-ID 030536. 



Most of our sources are selected from subsamples which are in 
principle not biased in favour of AGN. However, selecting the 
sample by using the availability of XMM-Newton data probably 
introduces a selection effect in favour of the presence of an AGN. 
Moreover, estimating the completeness level of this sample is 
difficult, sinc e it is not flux limited. F rom IR and sub-mm unbi- 
ased survevs. lRowan-RobinsonI (l2000h estimates that the number 
of HLIRGs brighter than 200 mJy at 60 fim over the whole sky 
is 100-200. 14 of them are included in our sample, which is the 
largest sample of HLIRGs studied in X-rays. 

Table [T] describes our sample. Column 2 shows the opti- 
cal spectral classification as derived from the literature: twelve 
sources in our sample present AGN characteristics. Eight of 
them are classified as 'type 1', and four of them as 'type 2'. We 
have classified QSO instead of Seyfert to those objects with in- 
trinsic 2-10 keV luminosity > lO^'* erg s ' (see Table|4]i. A cou- 
ple of sources have been classified from the literature as "Narrow 
line" (NL) sources, i.e. which show narrow forbidden emission 
lines typical of HII regions. All 'type 2' and one NL-SB galaxy 
are CT candidates. 

The analysis of the IR Spectral Energy Distribution (SED) of 
our sources has revealed that the SED can be modeled b y a com- 
bination of an AGN and a SB components (iRowan-RobinsonI 
I2000t IVerma et al.ll200l iFarrah et aljr2002ah . In Table [l] col- 
umn 3 we report the fraction of the AGN and SB component 
needed to fit the SED. Three objects are completely dominated 
by the AGN component, while in other three the SB component 
is dominant. 

In order to compare the properties of HLIRGs with other 
similar classes of objects, we have included in our study two 
samples. We cho se a s ample of 10 ULIRGs studied in X-rays by 
XMM-Newton dFranceschini et al.l l2003b . The sample is flux- 
limited at 60 i^m and complete to Seo nm > 5.4 J y. In addition, 
we ha ve selected all the HLIRGs (six) from the IStevens et al.l 
(1200 5*) sample of high redshift X-ray Compton-thin absorbed 
QSO. These sour ces have been observed in X-rays by RO SAT 
jPage et al.ll2"00l and by XMM- Newton dPage et al.ll2007l) . and 
by SCUBA in the sub-mm band (iPage et al.ll2004 . 



3. X-ray data reduction and analysis 

3.1. Data reduction 

Table|2]presents the most relevant information about the XMM- 
Newton observations. The data have been processed using the 
Science Analysis System (SAS) version 6.1.0, and have been 
analyzed using the standard software packages (FTOOLS and 
XSPEC) included in HEAsoft 5.3.1. 

We reprocessed the EPIC PN and MOS Observation Data 
Files (ODFs) to obtain new calibrated and concatenated event 
lists, using the SAS tasks EMPROC and EPPROC, including flie 
latest calibration files at the time of reprocessing. 

The new event files were filtered to avoid intervals of flar- 
ing particle background, and only events corre sponding to pat - 
tern 0-12 for MOS and 0-4 for PN were used tEhle et al.ll2005h . 
The events with energy above 12 keV and below 0.2 keV were 
also filtered out. The source spectra were extracted from circu- 
lar regions, whose radius was chosen in each case to optimize 
the signal-to-noise ratio (S/N), and to avoid the CCD gaps. The 
background spectra were taken in circular source-free regions 
near the object, also avoiding CCD gaps. We generated our own 
redistribution matrices and ancillary files (correction for the ef- 
fective area) using the SAS tasks RMFGEN and ARFGEN. 

XMM-Newton has detected 10 out of 14 sources (~ 70%) 
with different S/N quality. In those cases where the S/N ratio was 
poor, the MOS and PN spectra were co-added (Page et al. 2003|. 
All spectra were rebinned to have > 25 counts per energy chan- 
nel, except IRAS 00182-7112 (> 15 counts) and PG 1206+459 
(> 20 counts). The resulting EPIC spectra (see Fig.[Tli reveal het- 
erogeneous spectral properties for these objects (see Table|4]and 
Sects. I3.4l3.5l 

3.2. Non-detected sources 

We have estimated upper limits to the luminosity of the sources 
not detected by XMM-Newton. We estimated the count rate 
which would correspond to 3cr fluctuations of the background in 
a circular region of the PN-EPIC images, centered in the source 
coordinates. To convert between count rate and physical units a 
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Energy (keV) 

(c) PG 1206+459: PN, MOSl, MOS2 




Energy (keV) 

(d) PG 1247+267: PN, MOSl, MOS2 




Energy (keV) 

(e) IRAS F12509+3122: PN, MOSl, MOS2 




10-s 

0.5 


-0.5 



0.5 1 2 5 

Energy (keV) 

(f) IRAS 12514+1027: PN+MOS1+MOS2 



Fig. 1. XMM-Newton X-ray spectra and residuals of the detected sources from the HLIRGs sample. The solid Une is the best fit 
model. Continues on next figure. 
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Fig. 1. Continue 



simple model was chosen: a power law with F - 2 and Galactic 
absorption. 

3.3. Cluster emission subtraction 

Two sources of our sample, IRAS 09104+4109 and IRAS 
18216+6418, reside in clusters. They present soft extended 
emission from the intra-cluster medium (ICM). To take into ac- 
count this residual foreground in the subsequent spectral anal- 
ysis, we added an XSPEC thermal component to the spectral 
model. Two parameters characterize this model: temperature and 
normalization. 

We estimated the temperature of the cluster extracting an X- 
ray spectrum in an annular region around the source and fitting 
it with a thermal bremsstrahlung model. The normalization was 
obtained re-normalizing the flux from the annulus to that from 
the source circular region. To this end, we integrated the X-ray 
surface brightness profile of the cluster over both region^ 



Table 3. Chandra observations description. 



Source 


Instrument 


Grating 


Exp. time, 
[ks] 


Obs. date 


IRAS 09104+4109 


ACIS-S 


None 


9.17 


1999-11-03 


IRAS 18216+6418 


ACIS-S 


LETG 


171.82 


2001-01-18 



We determined the brightness profile using public Chandra 
data (see Table|3]l. Assuming an isothermal ICM, the radial pro- 
file of a cluster ca n be well fitted by a j6-model (Sarazin 1981 
Mushotzi^Hool). This profile was then convolved with the 
XMM-Newton PSF, before integrating it over the regions of in- 
terest. 



' The background spectrum of the annular region was extracted in a the XMM-Newton "blank fields" jRead & Ponmarill2003h to extract the 
circular region free of sources away from the cluster emission. We used background spectrum of the source. 
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3.4. Spectral analysis 

Our aim is to estimate the AGN and SB contribution to the total 
X-ray emission. The general model of the X-ray spectrum emit- 
ted by an AGN has typically four components: an underlying 
absorbed power lawQ, a reflection component, an iron K„ emis- 
sion line and a soft excess above the power law at energies below 
~ 1 keV. 

The typical power law photon index for AGN is T ^ 1.5 -2.1 

jNandra & Pounds 1994; Reeves & Turner 2000; Mateos et all 
I2OO5I). The soft excess c omponent, common in type 1 AGN 
dReeves & Turned 120001; iPiconcelU etall |2005|) . can be fit- 
ted with a thermal model, with an ex pected temperature re- 
markably con stant around 0.1-0.3 keV (' Piconcelli et al.l 120051: 
iGierlinski & Do ne 2006). The ratio between the luminosity of 
the soft excess component and the power la w component in 
flie so ft band (0.5-2 keV) is ~ 0.15 - 1.45 (iPiconcelh et all 
l2005h . The X-ray-to-bolometric luminosity ratio in the hard 
band (2-10 keV) for an AGN is ~ 0.015 (lElvis et al.l[T99l: 
iRisaliti & Eivisll2004 see Sectgli. 

Several models have been proposed for the origin of the ob- 
served soft excess in AGN, such as a relativistically blurr ed pho- 
toion ized disc reflection (Ross & Fabian 1993; Crumm y et al.l 
I2OO6I) or ionized absorption in a wind from the inner disc 
(iGierUiiski & Donell2004il2006h . It has also been shown that the 
analysis of the high resolution XMM-Newton RGS spectra is 
important to understand the nature of the soft excess emission 
(Bianchi et al. 2006). We have only used thermal models in this 
paper for simplicity, since more detailed models are not war- 
ranted by the quality and low spectral resolution of the data. A 
deeper analysis of the soft excess emission of these sources will 
be presented elsewhere. 

The SB emission can be modeled with a power law with 
a typical photon index F ^ 1-0 ~. 1-4 jWhite et all 119831; 
iDahlem et all 119981; iPersic & Rephaelil |2002|) . or with a ther- 
mal model whit tempe rature 0.5 - 1.0 keV (llwasaw il [19991 
iFranceschini et 31.12003 *). The soft X-ray-t o-bolome tric luminos- 
ity ratio for a SB is ~ 10 as found by Iwasawal ( 11999,) for a 
sample of four prototype powerful FIR SB galaxies. We can use 
this to estimate the bolometric luminosity of the sources. 

All sources were analyzed using the following scheme. First 
we fitted the data with a power law model plus Galactic ab- 
sorptiorQ (see TableH] column 3). We added intrinsic absorption 
(XSPEC zpha model) where it was statistically significanfl 

We compared this model with a power law reflected 
from neutral material (mod eled with the XSPEC pexra\|3, 
iMagdziarz & Zdziarskilll995h . When two models had the same 
number of parameters (and hence the F-test is not useful), the fit 
with the lowest was taken as our baseline model. However, 
when the from two models were comparable, we adopted as 
our baseline fit the one with less uncertainty in the determination 
of the model parameters. 

Then, we checked for any significant additional component: 
iron emission line (modeled with zgauss at ~ 6.4 keV) and/or 



soft excess. We fitted the latter with different thermal models 



^ '£ch where A' are the number of counts and e is the energy of 

the counts. 

^ All models referred in this paper include a multiplicative 
Galactic absorption comp onent fixed at the Galactic N/y value from 
iDickev & LockmaiJ ( fl990h . 

To this end we used the F-test, accepting additional spectral com- 
ponents only when they improved the fit with a significance > 3it. 

^ We kept fixed all the parameters of the pexrav model to the stan- 
dard values except the photon index of the incident power law, the re- 
flection scaling factor and the normalization. 



(blackbody: zbbody; bremsstrahlung: zbremss, Kellogg et al 
1975; br emsstrahlung with emission lines: mekal, Mewe et aL 
ll985i;lKaastra & MewellT993 ). We introduced this thermal com- 
ponent to parameterize the starburst emission and/or the AGN 
soft excess emission. 

In those cases where F was out of the range expected for 
AGN or SB, this parameter was fixed to 2, and we tried to fit 
again the spectrum with an absorbed power law and a reflection 
model (pexrav). In either case, additional intrinsic absorption 
and thermal emission models were also added (if needed) to im- 
prove the fit. 

The best fit model for each source is given in Table H] We 
have calculated the luminosity for each component in the hard 
and soft X-ray bands, coiTected by Galactic and intrinsic absorp- 
tion. A more detailed description of the analysis and the results 
for each source is presented in Section [33] 

We have estimated upper limits (see Table|4] column 11) for 
a thermal component in the X-ray spectra of those sources where 
it was not significant. To this end, we fixed all parameters to their 
best fit values. We added a thermal component (zbremss) with 
a fixed temperature, kT = 0.6 ke V (the mean temperature o f the 
ULIRGs thermal component from lFranceschini et al.ll20()3 l). and 
we calculated the 2cr confidence interval for the normalization 
parameter, which was then used to estimate the upper limit to 
the luminosity. 

3.5. Source by source analysis 
IRAS 00182-7112 

This type 2 QSO was detected only in the EPIC-PN camera. We 
modeled the spectrum of IRAS 00182 -7112 with a reflection 
comp onent, using the pexrav model jMagdziarz & Zdziarskil 
I1995 V The photon index is fixed to 2 and a pure reflection 
component is assumed. This implies a lower limit to the col- 
umn density of the absorber material (Nh > 10^^ cm"^). We 
marginally detect a narrow emission line at 6.75;^y j'j keV (sig- 
nificance 2 - 3cr), consistent with He-like Fe line (this energy it 
is also consistent at 2cr level with neutral Fe 6.4 keV line). 

ISO and Spitzer IR data sug gest the presence of a deeply 
obscured nuclear power source (iTran et al.l 1200 it I Spoon et aP 
2004). This is qualitatively consistent with our X-ray analysis 
results. The X-rays detected by XMM-Newton should be the re- 
flected emission from the AGN, and the line is an iron Kq, flu- 
orescent emission from the reflecting material. The equivalent 
width of this line, 0.8 ± 0.6 keV, is consistent with the CT hy- 
pothesis, but the poor quality of the spectrum prevents us from 
reaching any stronger conclusion. Assuming that the direct X- 
ray emission is completely absorbed by CT material, we have 
found that the intrinsic 2-10 keV luminosity of the AGN respon- 
sible for the reflection component is 6.3 x IQ'^^iln/Q.refi) erg s"' , 
where Q.refi < 27r is the solid angle subtended by the reflector at 
the illuminating source. This X-ray emission gives an estimate to 
the bolometric luminosity of the AGN of ~ 1.1 x lO'^L©, which 
is consistent with the IR observations (see Figs.|2]|3]l. 

In conclusion, our X-ray analysis of this source points to an 
AGN with CT obscuration. 



IRAS F00235+1024 

XMM-Newton observed this NL S B galaxy for 26 ks, b ut it was 
not detected by the EPIC cameras. IWibnan et al.l (l2003l) . assum- 
ing a thermal mekal model (kT = 0.5 keV), estimate an upper 
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Table 4. XMM-Newton spectral analysis results. Fluxes and luminosities in CGS units. Errors quoted are 90% confidence level for 
one interesting parameter throughout this paper. 



Source 


Model " 


Nh" 


r 


kT 


EW 










log 














(keV) 


(keV) 










J :>A 
TH 


J tiX 

'^PL 


J tlx 


IRAS 00182-7112 


D + E 


4.24 


2, E, = 9+f 




0.8 ±0.6 


5.0/9 


10.7 


-12.8 


44.9 


<41.9 


44.8 


<40.7 


IRAS F00235+I024 


- 


5.07 


2 


- 


- 


- 


- 


<-14.5 


<42.2 


- 


<42.4 


- 


IRAS 07380-2342 


- 


64.3 


2 




- 


- 


- 


<-13.7 


<41.7 


- 


<42.5 


- 


IRAS 09104+4109 


(B)' + A + E 


1.82 


1.62 ± 0.07 


- 


0.38!»:!? 
0.2if 


374/324 


42 


-11.7 


44.2 


<43.0 


44.5 


<41.8 


+BeppoSAX: 


(BY + B + D + E 


II 


'-■^-0.2 


T 1 +0.4 
-"•^-0.3 


330/326 


53 


II 


44.7 


44.2 


45.3 


44.2 


PG 1206+459 


A 


1.31 


1.7 ±0.9 






72/67 


- 


-12.5 


45.8 


<44.0 


45.1 


<42.8 


PG 1247+267 


A + C 


0.90 


1 98+005 

-097 






235/281 


30 


-12.3 


45.9 


45.5 


45.9 


44.0 


IRASFI2509+3122 


A + B 


1.24 




0.21 ±0.03 




38/30 


80 


-12.9 


43.8 


43.8 


44.3 


40.2 


IRAS 12514+1027 


B + F X Af 


1.67 


2 


35+" " 




5/15 


15 


-14.2 


43.2 


42.2 


43.3 


39.7 


IRAS 13279+3401 




0.99 


2 








<-14.2 


<42.1 




<42.2 




IRAS 14026+4341 




1.19 


2 










<-13.7 


<42.5 




<42.6 




IRASF14218+3845 


A 


0.93 


2.24 ±0.12 






70/75 




-13.1 


44.7 


<43.8 


44.6 


<42.5 


IRAS FI5307+3252 


A 


2.03 


2.1 ±0.4 






16/22 




-13.6 


45.4« 


<43.1 


45. 5« 


<41.9 


IRAS 16347+7037 


C + A 


4.48 


1.77 ±0.13 


1.53 ±0.18 




561/547 


26 


-11.7 


45.8 


45.7 


46.0 


45.4 


IRAS 18216+6418 


(C)' + C + A 


4.04 


1 C-7+0.10 


49+0.09 




329/332 


151 


-10.7 


45.2 


45.1 


45.6 


43.6 



" XSPEC models: A: power law, B: mekal, C: zbremss, D: pexrav, E: zgauss (emission line), F: zpha (intrinsic absorption). 

* Galactic neutral hydrogen column density (in units of 10"" atoms cm"^), from lDickev & LockmanI | |199C1|) . 
Observed frame 0.5-10 keV flux. 

'' Rest frame intrinsic X-ray luminosity. 
SX: Soft X-ray (0.5-2 keV) band; HX: Hard X-ray (2-10 keV) band; PL: Power law component; TH: Thermal (soft excess) component. 
'' This component models the cluster contribution. 
f Intrinsic Nh = 4+f x 10^^ cm^-^ 

* Luminosity corrected assuming iPanessa et al.l i2006l) results. See Sect. l3.5l for details. 



limit to the 0.5-2 keV luminosity of 2.8 x 10^^ erg s"', consistent 
with our result in this band. The limit in the soft band implies 
a weak SB emission in X-rays, compared to that expected from 
the IR luminosity (iFarrah et alj|2002al see Table |5j. Their upper 
limit to the hard X-ray luminosity (~ 1.9 x 10"*"^ erg s"') is also 
consistent with ours. 



IRAS 07380-2342 

This NL object has not been detected by XMM-Newton. 
IRAS 09104-1-4109 

The typ e 2 QSO IRAS 09104-f4109 resides in a rich cluster 
(Kleinma nnetan ri988). The X -ray soft ext ended emission from 
the I CM was already detected by ROSAT (iFabian & Crawford! 
Il995i) . We subtracted this foreground as explained in Section [33] 

The source spectra were extracted from a circular region of 
20" for all detectors, and the cluster spectra from an annular 
region between 20" and 90" (constrained to the CCD where 
the source is located). The cluster emission was fitted with a 
mekal model. The temperature is kT = 5.5 + 0.4 keV, and the 
metal abundance is 0.30 + 0.12 Z©, which is consistent with the 
mean Fe abundance for clusters with a temperature greater than 
5 keV dBaumg artner e t al]|200ll) . We obtained the radial bright- 
ness profile using Chandra data: a y6-model with a core radius of 
4" .6 + 0".6 and/? = 0".576!|;" °|^. 

The cluster emission represents 62% of the total 0.5-10 keV 
luminosity. The spectrum of the source was fitted with a power 
law with a Fe Ka broad (cr = 0.27 + 0.09 keV) emission line 
at 6.61;^Q jQ keV. This broad line could be explained as a com- 
plex of neutral and ionized narrow lines merged due to the low 
resolution of the detector. 



Although IRAS 09104-H4109 is classified as a type 2 QSO, 
the XMM-Newton data did not reveal any intrinsic absorp- 
tion feature. However, the Chandra observation of this source 
suggested a column density of 3 x 10^^ cm"^ (ilwasawa et al.l 
2001) . Also, BeppoSAX detected this object at energies greater 
than 10 keV, pointing to non-th ermal quasar emission e merg- 
ing from a thick absorbing torus jFranceschini et alll200d) . with 
Nh ~7x 10^'^ cm-2 . 

Our combined analysis of the BeppoSAX and XMM- 
Newton data sets (which will be taken as the best fit for this 
source in what follows) shows that a reflection-only model is 
needed to explain the complete spectrum in the 0.2 to 50 keV 
range. This implies a lower limit (Nh > 10^^ cm"^) to the 
column densit y of the absorber, which it is consistent with 
Ilwasawa et al.l (1200 ll) . where a similar analysis of this source 
is done with BeppoSAX and Chandra observations. Assuming 
r = 1.4, which is i n the flatter side of t he photon index dis- 
tribution of quasars ('Reeves & Turner ''2000), they found that a 
cold reflection model without transmitted component fits well 
the complete data. 

This model gives an intrinsic hard X-ray luminosity of 2 x 
lQ'^^(2n/Q.refd ergs"' for the AGN. The estimated bolometric 
luminosity of this source is then ~ 3.5 x lO'^L©, which is con- 
sistent with its IR data (see Figs.|2][3]l. 

We have also found in our combined analysis a new thermal 
component, but its temperature (kT ~ 3 keV) is too high to be 
associated to SB or AGN emission. Since ther e is no evidence of 
a SB component in the IR SED of this source (i Rowan-RobinsonI 
2000), this thermal component is probably due to an incomplete 
subtraction of the cluster emission. Previous results indicate that 
a strong cooling flow of the ICM is taki ng place in the core 
of the clu ster (Fabian & Crawford 1995; Allen & Fabianlll998t 
lEttori & Fabian. ,1999; .Iwasawa et al.. ,20011) . so the isothermal 
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ICM hypothesis we have assumed could underestimate the clus- 
ter luminosity in the central region. 

We can confirm that IRAS 09104+4109 is probably a CT 
source, and that the X-ray emission detected below 10 keV by 
XMM-Newton is only a reflection continuum from cold matter 

PG 1206+459 

The XMM-Newton observation of this QSO is contaminated by 
background flares at the beginning and at the end of the obser- 
vation. The final effective exposure is ~ 7 ks. The spectra con- 
tinuum has been modeled with a power law. We have not de- 
tected significant intrinsic absorption or soft excess. The lack of 
absorption in the X-ray spectra is consistent with th e optical and 
IR evidences (iHaas et al.ll99 8; Rowa n-Robinso n 2000). Its hai'd 
X-ray luminosity is 1.3 x 10^^ erg s"', which gives a bolomet- 
ric luminosity of ~ 2.2 x IO'-'Lq, in agreement with the IR data 
(see Figs. |2l [3]). The X-ray spectrum of this object is therefore 
consistent with having a pure AGN origin. 

PG 1247+267 

The X-ray spectrum of this QSO has been modeled as a power 
law and a thermal component. A pexrav model is formally the 
best fit (x^/v - 224/282) of these data. However, the photon 
index obtained (F ~ 2.3) with this model is slightly larger than 
the one expected for an AGN. Moreover, the reflection scaling 
factor (~ 4) is quite larger than the typically expected for type 1 
sources (within and 1). No other reflection features have been 
found in the X-ray spectrum. Therefore, we have adopted the 
power law plus thermal component as our best fit. The bolomet- 
ric luminosity of the AGN is ~ 1 .4x IO'^^Lq, using its hard X-ray 
luminosity. This result is consistent with the IR observations (see 
Figs. E IS. 

The temperature of the thermal component {kT - 
0.48;^Q j^ keV) is consistent with the typical temperature of a SB 
galaxy. However, the bolometric luminosity that we can estimate 
through the soft X-r ay emission for th i s SB is ~ 10^^ erg s 
much higher than the lRowan-RobinsonI (1200 0*) estimate (< 5.2 x 
10^^ erg s '). Therefore, the soft excess component is too lumi- 
nous to have a pure SB origin. Furthermore, its soft excess-to- 
power law soft X-ray luminosity ratio is ~ 0.4, which is typical 
for soft excess observed in AGN. The X-ray spectrum of this 
source is consistent with being dominated by an AGN. 

IRAS F1 2509+31 22 

A significant fraction (~ 50%) of the observation time of this 
QSO is affected by high background. The PN and MOS spec- 
tra can be fitted by a power law model and a thermal compo- 
nent with kT = 0.21 + 0.03 keV, at a lower energy than that 
expected for a standard SB, but consistent with soft excess orig- 
inated in an AGN. The bolometric luminosity inferred from the 
X-ray thermal luminosity is one order of magnitude greate r than 
the estimate for a SB using IR data (iFarrah et al.ll2002ah . The 
thermal-to-power law luminosity ratio is ~ 1.2, in the range 
of AGN soft excess. The X-ray spectrum of this source is also 
AGN-dominated. 

IRAS 12514+1027 

This Seyfert 2 galaxy was also observed by ROSAT 
dWiknan et all [T998i) . but only XMM-Newton has been able 



to detect its X-ray emission. IWilman et alj (l2003h considered 
only the PN spectrum, modeled with a reflection component 
(pexrav, with fixed F = 2) and a thermal component (mekal, 
kT - 0.31;^Q Qg keV) corrected by intrinsic absorption (Nh = 
1.3+°-^ X lO^' cm-^). The resulting vis 8.5/12. 

In our analysis, the three EPIC spectra were coadded (see 
Sect. 13.11 ). We modeled the spectrum with an absorbed power 
law with a photon index fixed to 2, and a thermal component. 

The temperature of the thermal component is kT = 
0.35^QQ^ keV, consistent with that usually observed in a SB 
galaxy. The soft component-to-power law luminosity ratio is 
slightly lower than that associated to an AGN. The 0.5-2.0 keV 
luminosity of this soft component is 1.5x10'*" erg s"', which im- 
plies a bolometric luminosity for the SB of ~ 1 .5 x 10^^ erg s"' . 
The IR luminosity of the SB estimated from th e analysis of the 
IR SE P of the source is ~ 5 x 10"*^ erg s ' CRowan-RobinsonI 
l2000h . This thermal component could be consistent with a SB 
origin. 

The optical, IR and X-ray data point to a CT quasar and a SB 
of comparable bolometric luminosities. 

IRAS 13279+3401 

This QSO has not been detected by XMM-Newton. 
IRAS 14026+4341 

The XMM-Newton observation of this QSO 1.5 is heavily con- 
taminated by background flares. All the PN data are affected by 
count rate background greater than 15 counts per second. The 
MOS data have a brief interval free of flares, but the source is 
not detected. 

IRAS F1 421 8+3845 

The QSO IRAS F14218+3845 was observed by XMM-Newton 
in two occasions. The second observation was heavily affected 
by high radiation background. We co-added the six spectra of the 
different observations to increase the S/N ratio. The spectrum 
was modeled with a power law. No significant soft excess or 
intrinsic absorption was found. 

The IR data suggest that this HLIRG is a SB dominated 
source dVerma et alj|2002l: iFarrahet al.|[2002ah . but our analy- 
sis of its XMM-Newton X-ray spectrum does not reveal any 
SB features. Using the upper limit that we have estimated (see 
TableHli for a thermal component, the total SB luminosity is less 
than 6 x 10^^ erg s~', which i s consis tent with the SB luminos- 
ity estimated by Farrah et akl (l2002al) through IR and sub-mm 
data (6 x 10'*^ erg s~'). Although a SB component cannot be 
excluded, the data point to an AGN-dominated X-ray emission 
dFranceschini et al.ll2003h . 

IRAS F1 5307+3252 

Previous observations with R OSAT and ASC A d etected no X- 
ray emission from this QSO 2 dFabian et al.' 1996; O gasaka et alJ 
1 1997.) . We detected a faint X-ray emission in the XMM-Newton 
public data. The observation of this source is affected by high 
background flares. The three EPIC extracted spectra were coad- 
ded to increase S/N ratio. We fitted this spectrum using a power 
law. We were not able to find any absorption feature or thermal 
emission. 
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Table 5. Flux and luminosity data of the sample. 









Flux density (Jy]^ 




IR 




lo^ 


'L (erg s 






Source 


Rl 


lljim. 


25yum 


60y[(m 


100//m 


ModeF 






LlRJiR 


tAGN 
IR RR 


IRRR 


IRAS 00182-71 12 


17.738 


<0.06025 


0.133±0.010 


1.20+0.08 


1.19+0.12 


S+A 


46.49 


<46.72 


<46.93 


<46.48 


46.74 


IRASF00235+1024 


>21.5 


<0.173 


<0.193 


0.43+0.06 


<0.94 


S+A 


<46.76 


<47.27 


46.74 


46.44 


46.45 


IRAS 07380-2342 


16.869 


0.48±0.03 


0.80+0.08 


1.17+0.09 


3.5+0.3 


A+S 


46.56 


47.08 


46.97 


46.79 


46.48 


IRAS 09104+4109 


17.819 


0.13±0.03 


0.334+0.013 


0.53+0.04 


<0.44 


A 


<46.42 


<46.99 


<46.92 


46.84 


<46.15 


PG 1206+459'' 


15.135 


0.21±0.04 


<0.113 


0.26+0.05 


0.35+0.07 


A 


47.20 


<47.94 


<47.80 


47.78 


<46.57 


PG 1247+267 


14.621 


<0.126 


<0.113 


0.24+0.05 


0.17+0.03 


A 


47.70 


<48.39 


<47.94 


47.91 


<46.76 


IRASF12509+3122 


16.590 


<0.106 


0.10+0.03 


0.22+0.04 


<0.675 


A+S 


<46.86 


<47.37 


47.00 


46.76 


46.62 


IRAS 12514+1027 


17.654 


<0.0632 


0.190+0.016 


0.71+0.06 


0.76+0.15 


S+A 


46.18 


<46.51 


46.63 


46.27 


46.39 


IRAS 13279+3401 


15.689 


<0.0937 


<0.126 


1.18+0.08 


1.20+0.18 


A+S 


46.58 


<46.84 


46.53 


46.36 


46.04 


IRAS 14026+4341 


15.651 


0.12±0.03 


0.285+0.014 


0.62+0.06 


0.99+0.24 


A+S 


46.26 


46.70 


46.54 


46.34 


46.11 


IRASF14218+3845 


>21.5 


<0.0969 


<0.075 


0.57+0.06 


2.10+0.17 


S+A 


47.80 


<48.06 


46.86 


46.15 


46.76 


IRAS F15307+3252 


19.131 


<0.065 


0.071+0.019 


0.23+0.04 


<0.71 


A+S 


<47.07 


<47.46 


47.22 


47.05 


46.73 


IRAS 16347+7037 


13.979 


0.059+0.010 


0.122+0.004 


0.27+0.05 


0.35+0.07 


A+S 


47.42 


47.86 


47.81 


47.73 


47.04 


IRAS 18216+6418 


13.403 


<0.238 


0.445+0.012 


1.24+0.05 


2.13+0.17 


A+S 


46.49 


<46.89 


46.78 


46.54 


46.37 



" Observed by IRAS (from NED). 

* UK-R magnitude (SuperCOSMOS Sky Survey). 

AGN (A) and/or starburst (S) components needed to fit the IR SED (as in Table[T] col. 3). First letter indicates the dominant component. 

Infrared luminosities in the 40 - 500 fiin (FIR) and 1 - 1000 jiin (IR) bands, computed using IRAS fluxes dSanders & Mirabellll996l) . 



IR dHaas et al. l ll998HFarrahetalJl2002bl) observations suggest. 
Previous X-ray data from ASCA was also consistent, and there 
was no evidence of iron Ka emission feature or any absorption 
edge dNandraetal .111995b . 

The soft excess has L(0.5 - 2.0 keV) = 5.6 x lO'*^ erg s"'. 
This would imply a SB bolometric luminosity three orders of 
magnitude greater than the SB luminosity calculated with the IR 
data (Farrah et al. 2002b, see Table|5]l. Therefore this component 
is unlikely to be associated to a SB. Moreover, its thermal-to- 
power law luminosity ratio is consistent with a soft excess from 
the AGN. 

This model gives a bolometric luminosity for the AGN of 
~ 2.2 X IO'^Lq, in agreement with the IR observations of this 
source (see Figs.|2]|3]l. Our analysis points to a pure AGN origin 
for the X-ray spectrum, in agreement with the optical and IR 
data. 



IRAS 18216+6418 

The PN data of this QSO 1 .2 were heavily affected by pile-up, so 
we used only the MOS data. The source spectra were extracted 
from a 20" radius circular region in both MOS detectors. This 
QSO is located in a rich cluster, and ROSAT detected the ICM 
thermal emission (Saxton et al. 1997; Hall et al. 1997). We sub- 
tracted the soft X-ray emission from the cluster, as explained in 
Sect. [33] 

The MOS cameras operated in small-window mode, so we 
considered the PN image to model the cluster (the pile-up only 
affects the central region of the source). We extracted a spectrum 
from an annular region between 20" and 80" (constrained to the 
CCD where the source is located). The resulting temperature of 
the zbremss model was kT = 2.3!i;° keV We used also th e 
radial X-ray brightness profile published by iFang et aTl ( |2002|) . 
(core radius of 17".6 + 0".17,/3 = 0".74+o" °^), to renormalize 
the cluster model. The cluster X-ray emission is 32% of the total 
0.5-10 keV luminosity. 

The source spectrum best fit is a power law (F - 1.57^y jj) 
with a soft thermal component. A pexrav model is formally the 
best fit ix^/v = 321/333) of this spectrum. As discussed in PG 



60 and 100 fim fluxes are ISO data from iHaas et ah t200Q) . 



XMM-JVewtoJi has observed IRAS F15307 +3252 on two 
more occasions, but the data are still private, 'iwasa wa et alJ 
d2005h . using the complete data set, found a prominent Fe Ka 
line at ~ 6.5 keV, indicating the presence of a CT AGN. This 
is in agreement with optical spectropolarim etry data indicatin g 
the presence of a dust-enshrouded quasar dHines et al.lfT995h . 
The estimate of the AGN bolometric luminosity using the ob- 
served emission line luminosit y is also cons i stent with previous 
result s lYun & ScovilldT99l lAussel et alJ Il998l: IVerma et alJ 



ts r 

I2OO2I; [Peeters et al.i.2004l) . The hard X -ray emission detected by 
us is probably reflec ted radiation, because of CT obscuration. 
iPanessa et al.l d2006h found that the ratio between the intrinsic 
and the observed X-ray luminosity in CT Seyfert galaxies is 
~ 60. We have corrected the X-ray luminosity calculated with 
the public XMM-Newton data by this factor The resulting hard 
X-ray luminosity (~ 3.2 x 10"^^ erg s~ ') is consistent with the es- 
timate given by I wasawa et al.l d2005') (L^ > 1 x lO**^ erg s"' ), us- 
ing the luminosity of the iron emission line. Assuming an AGN 
origin, this X-ray emission gives a bolometric luminosity for this 
source of ~ 5.5 x lO'^L©, in agreement with the IR data (see 
Fig s.lZlllll. 

llwasawa et al" I d2005h also found extended soft emission, 
with kT - 2.1;^Q4 keV. They identify this extended emission 
with hot gas associated with a relatively poor cluster around this 
object. Although no galaxy cluster has been found associated 
to this source, an HST observation shows a moderate galaxy 
over-density. Moreover, its bolometric luminosity to tempera- 
ture relation would be similar to that typical of poor clusters 
dFukazawa et al.ll2004 . 

In summary, the X-ray spectrum of this source is consistent 
with the emission originated in a CT AGN. 



IRAS 16347+7037 

The spectrum of flie QSO IRAS 16347 H-7037 was modeled with 
a power law and a thermal component. No intrinsic absorption 
is detected. 

The XMM-Newton X-ray spectru m is consistent w ith a 
type 1 AGN spectrum, as the optical (lEvans et alJ Il998l) and 
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1247+267, the steeper photon index (F ~ 2.3) and the reflection 
scaling factor » 1 (7? ~ 15) lead us to adopt the power law plus 
thermal component as our best fit. 

The photon index of the power law is not consistent with 
previous X-ray obser vations with ASCA (F = 1. 75 ± 0.03, 
Yamashitaetal]ll997h and Chandra (F = 1.761+g ;jg. lFang et al.1 
2002 V 

In Fig |l(j)| a systematic effect in the Ax^ spectrum can be 
seen above 5 keV. This may be due to a bad extraction of the 
cluster emission or to the use of the blank-field background. If 
we ignore the data above 4.5 ke V, we get a steeper photon index 
(F - 1.68 + 0. 11), compatible with previous X-ray observations. 

A thermal component is also detected, with a temperature 
of 0.49;^y IJg keV, consistent with SB emission. However, if this 
emission were associated to the SB, the bolometric luminosity 
of the SB would be three orders of ma gnitude higher than t he lu- 
minosity calculated using the IR data jFarrah et alj|2002al) . The 
soft component-to-power law luminosity ratio is in the range of 
that typicall y observed in an AGN. 

Ginga dKii et al.' '1991), ASCA ("Yamashita et al.' '19971) 
and Chandra (Fang et al. 2002; Yaqoob & Serlemitsos 
2005h detected iron emission features in this HLIRG. 
Jimenez-Bailon et a D (I2007i) detected a Fe-K emission line 
with a complex structure in the FN XMM-Newton spectrum 
of this source. We have also detected an emission line in the 
6-7 keV rest frame band, but the significance of the detection 
(< 2cr) was below our adopted threshold, and therefore we 
have not considered it further. We have estimated a 3cr flux 
upper limit to a broad (cr = 0.1 keV) line component at 6.4 keV 
of < 3 X 10"^ photons cm"^ s"\ consistent with the value 
~ (3 + 1) X 10"^ p hotons cm"^ s ' obtained with the Chandra 
data (iFang et al.ll20 02). 

Assuming an AGN origin, the hard X-ray emission gives a 
bolometric luminosity for this source of ~ 7 x IO'-'L©, in agree- 
ment with its IR data (see Figs.|2][3]). The X-ray spectrum of this 
object is consistent with a pure AGN origin. 

4. Discussion 

In their study of X -ray emission from ULIRGs, 
iFranceschini et al.] (l2003h define that the X-ray emission 
of a ULIRG is AGN-dominated if it presents either: a) a high 
X-ray luminosity, L(2 - 10 keV) > 10"*^ erg s"'; b) a heavily 
obscured hard X-ray component with Nh > 10^^ cm"^ (very flat 
or inverted hard X-ray spectra); or c) a Fe-K emission complex 
at ~ 6.4 keV with equivalent width > 1 keV (iron fluorescent 
emission from material illuminated by the AGN). The ten 
detected HLIRGs from our sample present at least one of the 
three characteristics above, thus showing an AGN-dominated 
X-ray spectrum. 

This resu l t is i n agreement with the trend noted by 
IVeilleux et all ([T999I) for ULIRGs: the fraction of sources with 
Seyfert characteristics increase with L/s (from ~ 25% among 
ULIRGs with L/r < IO'^-^Lq to ~ 50% among those with 
L/s > 10'^ "^Lq). However, we must keep in mind that our sam- 
ple is not complete, and could be slightly biased towards AGN. 
The first subsample in Rowan-Robinson (2000), which is not bi- 
ased in favour of AGN (see Sect. |2]l, presents 50% of objects 
with Seyfert characteristics. 

Five objects from our sample (four type 2 objects and one 
NL object) are probably CT, as reported in the literature. Our 
analysis of the XMM-Newton data of two sources (IRAS 00182- 
7112, IRAS 12514+1027) are consistent with the CT hypoth- 
esis, as well as our combined analysis of XMM-Newton and 
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Fig. 2. 2-10 keV X-ray luminosity of the power law compo- 
nent versus FIR luminosity (40-500 /vm). Filled symbols repre- 
sent sources where we have detected X-ray absorption (models 
D or F in Table HJ. The top thick solid line and area between 
the thin solid lines ("AGN-zone") indicates the X-ray luminos- 
ity e xpected for a n AGN of a given FIR luminosity (90% disper- 
sion, ElvEiliDIilM; jiiMlM-&EMs 2^^ The bottom thick 
solid line ("SB-line") indicates the X-ray luminosity expected 
for a SB of a given FIR luminosity dPersic et alJl2004l ; lKennicut3 
Il998h . See Sect. |4]for details. 



BeppoSAX data for IRAS 09104+4109. IRAS F00235+1024 
has not been detected, probably because it is heavily obscured. 
We have found no absorption features in IRAS F15307+3252, 
but as explained in Sect. 13.51 recently published results from 
XMM-Newton private data are consistent with CT absorption 
(Iwasawa et al. 2005). 

We have calculated th e FIR (40-500 /zm) lumino sities (L/t/r) 
using the IRAS fluxes (ISanders & Mirabell 1 19961 see Table 
|5] column 7). In Fig. |2] we have plotted the 2-10 keV lu- 
minosity of the power law component for each source ver- 
sus the FIR lumi nosity. We ha ve included the ULIRGs data 
from IFra nceschin i et al.l (l2003b and the high-z QSO from 
IStevens et al.i (i2005 ^ for comparison. No correlation between the 
2-10 keV and the FIR luminosity is found in HLIRGs, although 
it must be kept in mind that the sample is not complete in any 
sense. 

We have estimated the expected X-ray luminosity for a stan- 
dard AGN, given its FIR luminosity. To this end, we calculated 
the 2-10 keV-to-FIR luminosity ratio typical of ne arby bright 
QSOs, using their mean SED. Since the original lElvis et alj 
( 1994) sample was biased tow ards sources with high X-ray lumi- 
nosity, we have employed the iRisahti & Elvis! (2004) new data 
on QSO SED. Using this corrected SED, the ratio of the 2- 
10 keV band to bolometric luminosity changes significantly by 
a factor of 2, from ~ 0.03 to ~ 0.015. The X-ray luminosity de- 
rived from the IR luminosity using the latter ratio is plotted in 
Fig. |2] with a thick solid line. The top area between thin lines 
("AGN-zone") is the dispers ion of the SED , calculated with the 
90 percentile distribution (Elvis et al.lll994 . 

We can also calculate a relationship between FIR and X- 
ray luminosity for SB galaxies. The SFR of a SB is given by 
its F IR luminosi ty by SFRf/^ ~ i™/2.2 x 10"*^ Mq yr^' 
(Ke nnicuttlll998l) . and by its 2-10 keV X-ray luminosity by 
SFRx ~ L2-io/10^^ Mo yr~' (iPersic et alj|2004 afl luminosities 
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Fig. 3. AGN to total- IR luminosity ratio, calculated by 
iRowan-Robinsod (12000 '). versus AGN (estimated through our 
X-ray data) to total (computed with the IRAS fluxes) IR lu- 
minosity ratio. The dotted line means equal values for the ra- 
tios. In case of upper limits we plot as er ror bars the disper- 
sion of X-ray-to-IR luminosity ratio from lElvis et aL I (119941) . 
IRAS F00235+1024, IRAS 07380-2342, IRAS F12509+3122, 
IRAS 125 14-H 1027, IRAS 132794-3401, IRAS 140 26-H4341 and 
IRAS F14218-H33845 IR data are taken from the iFarrah et"!!] 
(1200 2a') results; IRAS F15307+3252 IR data is taken from the 
IVerma et al. (20 0j) results; IR AS 16347-1-7037 IR data is taken 
from the iFarrah et alj (l2002bl) results. Symbols as in Fig.|2l 



are in CGS units). Assuming equal SFR, the 2-10 keV to FIR 
luminosity ratio is L(2-10 keV)/Lf/R ~ 4.5 x 10"^. This relation 
is shown in Fig.|2]as the lower thick solid line ("SB-line"). 

Most HLIRGs and all high-z QSO are in the "AGN-zone", 
while only the AGN-dominated ULIRGs and two HLIRGs seem 
to be composite sources: their X-ray luminosity is too high to 
be produced only by a SB (above the "SB-line"), and their 
FIR luminosity is too high to be produced only by an AGN 
(to the right of the "AGN-zone"). The SB-dominated ULIRGs 
are concentrated near the "SB-line". The upper limits for the 
X-ray-undetected HLIRGs seem to indicate that their hard X- 
ray emission comes only from SB activity. However, the X- 
ray emission of the non-detected sources could be affected by 
heavy obscuration. Actually, one of the non-detected sources, 
IRAS F00235-H1024, is probably CT as seen from their IR 
spectrum so its X-ray emission could be depressed by heavy 
absorption. However the remaining sources are optical QSOs, 
where little or no absorption is expected. For example, the 
QSOs from IStevens et al.l ( l2005h show relatively low absorp- 
tion (21 < log Nh < 22 cm"^). Recent XMM-Newton observa- 



tions of these high-z HLIRG s suggest highly i onized winds with 
22.5 < logA^H < 23.5 cm'^ (iPage et al.ll2"007h . Furflier sensitive 
data on isotropic indicators (such as FIR or MIR or > 10 keV 
emission) are needed to investigate the seemingly contradictory 
nature of these HLIRGs. 

Note that if the lElvis et al.l (1 1 9941) X-r av-to-FIR ratio is 
used instead of the Risaliti & ElvisTi 20041) ratio, only three 
HLIRGs (IRAS 09104-1-4109, IRAS 16347-H7037 and IRAS 
18216-1-6418) would lie on the "AGN-zone", and the rest would 
be considered a s composite AGN/SB s ources. This confirms the 
relevance of the lRisahti & ElvisI (l2004i) correction. 



' F14218 FI5»' 



I 10 100 1000 
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Fig. 4. Lx (0.5-10 keV) of the soft excess component versus 
Lx (2-10 keV) of th e power law component. Asterisks repre- 
sent SB galaxies from lFranceschin i et al. (2003). The dotted line 
is a correlation obtained bv iFranceschini et alj (1200 3) for SB- 
dominated ULIRGs. The dash-dotted line is a correlation for the 
SB-dominated ULIRGs and the SB galaxies samples calculated 
by us. The soft excess in IRAS 09104-1-4109 is probably associ- 
ated to cluster emission (see Sect. |3.5| for details). Symbols as in 

Figm 



The HLIRGs FIR emission is systematically above the 
iRisaliti & ElvisI (|2004|) estimate for a standard local QSO of the 
same X-ray luminosity (i.e., the sources are located on the right 
of the thick upper line in Fig. |2]). This FIR excess could be as- 
sociated to the SB activity in HLIRGs. Alternatively, this could 
also hint to a possible difference between the standard QSO SED 
and the AGN component of the HLIRGs SED. In this line, it has 
been shown that the shape of the SED is probably related to the 
luminosity dMarconi et alj|2004l) . 

We tried to unravel the origin of the excess infrared emission 
with respect to that predicted using the Risaliti & EJvis| (|2004|) 
QSO SED. We have estimated the AGN contribution to the total 
IR luminosity of the HLIRGs in two independent ways, in order 
to know if this IR excess comes from SB activity, or from an 
intrinsic differenc e in the AGN SED . 

On one hand, Rowan-Robinso 3 (12000) has modeled the IR 
SED of all sources in this sample with radiative transfer mod- 
els, and he estimated the contribution of the AGN dust torus 
(^/« rr) *s t° *s luminosity (1-1000 yum, Ljr^rr) 

(Table|5] columns 9 and 10 respectively). We have corrected the 
relative conti'ibution of IRAS F00235-t-1024, IRAS 07380-2342, 
IRAS F12509-1-3122, IRAS 12514-1-1027, IRAS 13279-H3401, 
IRAS 14026+4341, IRAS F14218-1-3845, IRASF15307+3252 
and IRAS 16347+7037 with the results of Ver ma etall ( |2002|) . 
Farrahetal. (2002^) and iFarrah et al. (2002b), where new IR 
and sub-mm data are used. 

On the other hand, we have calculated this ratio using the 
power law X-ray luminosity. We assume tha t this emission is as - 
sociated to an AGN with a standard SED (Risaliti & Elvis 2004, 
and we calculate the IR-to-2-10 keV luminosity ratio, similarly 
to that done in Fig.|2] We obtain an independent estimate of the 
expected IR luminosity of the AGN (Lf^^^^). 

In Fig. |3]we compare the relative c ontribution of the AGN t o 
the total IR luminosity calculated by iRowan-RobinsonI ( l2000l) . 
to our estimate obtained from X-ray data (we have included only 
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Fig. 5. Lx (0.5-10 keV) of the soft excess component versus FIR 
luminosity. The soft excess in IRAS 09104-1-4109 is probably 
associated to cluster emission (see Sect. l3.5l for details). Symbols 
as in Fig|2] 

the X-ray detected sources). The dotted line corresponds to the 
1:1 relation, i.e. an agreement between the two estimates. Most 
of the sources have lower limits on the abscissae because their 
12 and 25 jjm IRAS fluxes are upper limits (see Table|5]l. 

Our estimates of the AGN relative contribution for all 
source s are formally consistent with that of Rowan-Robinson 
(I2000I) . However, there seems to be a systematically overesti- 
mate of the IR AGN component from iRowan-RobinsonI (|2000|) 
with respect to the X-ray measurements. 

The values of L^^^^^/Ljjf plotted in Fig. [3] are independent 
of the SB luminosity, so this disagreement is probably due to 
the IR-to-X-ray ratio used to estimate the IR luminosity from 
the X-ray luminosity. This favours the hypothesis of an intrinsic 
difference between the standard QSO SED and the AGN compo- 
nent of the HLIRGs SED. A detailed radio-to-X-ray study of the 
spectral energy distribution of HLIRGs is needed to solve this 
question. 

Fig. ID is a luminosity-luminosity plot of the soft X-ray 
excess versus power law components. In IRAS 09104-1-4109, 
IRAS F12509+3122, PG 1247 -h267, IRAS 16347+7037 and 
IRAS 18216-1-6418 the soft excess component is too luminous 
to come only from a SB (see Sect. l3.5l for a detailed description 
of each source). In the case of IRAS 09104-1-4109 the soft ex- 
cess is probably due to an incomplete subtraction of the cluster 
emission, while in the remaining four sources is probab ly of the 
same origin as in luminous QSO jPiconcelli et al.ll2005l) . 

An X-ray thermal component ass ociated to SB emis s ion ha s 
been observed in all ULIRGs from iFranceschini et al.] (|2003|) . 
However, only in 1 out of 14 HLIRGs we have found a soft X- 
ray emission whose origin could be associated to SB activity. 

Oddly, the above HLIRGs with AGN-like soft excess emis- 
sio n follow the correlation f ound for SB-dominated ULIRGs 
by IFranceschini et al.1 ( l2003l) (dotted line in Fig. |4]i. To in- 
crease the statistics and to test if this correlation holds at 
low er IR luminosities, we ha ve included a sample of SB galax- 
ies (IFranceschini et all l2003 l and references therein). We have 
calculated a non-parametric correlation coefficient (generalized 
Kendall's TaiQ) for the SB galaxies and the SB-dominated 

W e employed the ASURV software for this test jlsobe et al]| 19851 . 



ULIRGs, finding Zr = 3.69 with a significance leveQ of 99.98% 
(> 3cr). The correlation slope is consistent with that obtained by 
Franceschini et al. (2003) for SB-dominated ULIRGs only. 

We have investigated the relationship between the X-ray soft 
excess component luminosity and the FIR luminosity (Fig. |5]l, 
finding no clear correlation. A test using generalized Kendall's 
Tau confirms this impression: - 2.13, with a significance 
level 96.66% (< 3cr). 

We have also tested a possible cosmological evolution in the 
sample. We have estimated the SFR for each source using its 
IR luminosity (Kennicutt 1998). Their SFR and 2-10 keV-to-IR- 
luminosity ratio versus redshift are plotted in Fig.|6] Cosmic star 
formation sh ows an important decline between z ~ 2 and the 
present day (IFranceschini et all 1 19991) . so we expect an incre- 
ment of the SFR of HLIRGs up to z ~ 2. Higher SFR at higher 
redshift is observed in the upper panel of the Fig. |6] However 
the sources follow the lower envelope, which is the IRAS FSC 
sensitivity limit (solid line in Fig.|6]l, indicating clearly a selec- 
tion effect. Therefore, we can not draw conclusions about the 
dependence of SFR with redshift. 

As shown in the bottom panel of the Fig. |6] the ratio of hard 
X-ray-to-FIR luminosity remains constant with z- This holds 
even if we subtract the FIR luminosity emitted by the AGN, cal- 
culated using the X-ray PL luminosity as in Fig. [3] 

We have seen that the IR emission is consistent with an AGN 
origin, but if we assume that the IR excess shown in Fig.|2]is as- 
sociated to the SB activity, Fig.|6]shows that its evolution must be 
similar to that of the X-ray emission. This, in turn, suggests that 
the presence of a SB and the occurrence of AGN activity through 
accretion onto a super-massive black hole are physically related. 
This result is in agreement with the coeval black hole/stellar 
bulge formation hypothesis (iGranato et al.ll2004t IStevens et al.l 
120051: iDi Matteo et al.ll2005l) . 



5. Conclusions 

We have performed a systematic X-ray study of a sample of 14 
Hyper-Luminous Infrared Galaxies using XMM-Newton data 
from the archive, and our own private data. We modeled the X- 
ray spectra of each source, finding very heterogeneous spectral 
properties. Our results are summarized as follows: 

1 . All X-ray detected HLIRGs of the sample (ten sources) have 
AGN-dominated X-ray spectra. 

2. No correlation is found between the 2-10 keV and IR lu- 
minosities in HLIRGs. The hard X-ray luminosity of most 
(eight) HLIRGs is consistent with emission from an AGN 
component only. However two HLIRGs (as well as all 
the AGN-dominated ULIRGs) seem to show a composite 
AGN/SB nature: their X-ray luminosity is too high to be pro- 
duced only by a SB, and their IR luminosity is too high to be 
produced by only an AGN. The remaining four HLIRGs are 
undetected in X-rays. 

3. The hard X-ray luminosity associated to the AG N is system- 
atically below the one expected for a local QSO (Elv is et all 
.1994,: ,Risaliti & Elvis 2004 ) of the same IR luminosity. This 
seems to suggest that there is an intrinsic difference between 
the AGN component of the HLIRGs SED and the SED of 
local QSOs. A detailed radio-to-X-ray study of the HLIRGs 
SED is needed to understand this issue. 



' Note that even excluding the isolated source in the bottom left cor- 
ner, this correlation remains almost unchanged. 
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Fig. 6. Symbols as in F ig|2l Top: SFR derived from FIR lumi- 
nosity (iKennicuttll 1 998;b . versus redshift. The dotted line repre- 
sents the SFR limit corresponding to the HLIRGs definition. The 
solid line is the SFR for a sourc e with IRAS fluxes equals to 
the IRAS FSC sensitivity limits dMoshir et all 1 1990) . Bottom: 
2-10 keV power law-X-ray- to-FIR luminosities ratio versus red- 
shift. The area between the dotted lines represent the expected 
ratio for a quasar (with 90% of dispersi on) with a standard SED 
dElvis et al.lll994tlRisahti & Elvisll2004l) . 



4. There are five Compton-thick candidate sources in our sam- 
ple, including all the type 2 AGN (four) and one starburst ob- 
ject undetected in X-rays. We have found evidence for X-ray 
absorption in three of the type 2 AGN. Our analysis supports 
the Compton-thick nature of these sources. 

5. We have detected five HLIRGs with soft excess emission. 
Four sources present luminosities consist ent with type 1 
AGN soft excess. While all ULIRGs from Franceschi ni et alJ 
([2003) present a soft X-ray thermal component associated to 
a SB, only in one HLIRGs we have found a soft excess prob- 
ably associated to SB emission. 

6. The hard-X-ray-to-FIR luminosity ratio remains constant 
with z, suggesting that the AGN and SB phenomena are 
physically connected in HLIRGs. 
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